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Based on macroscopic quantum electrodynamics in linearly and causally responding media, we 
study the local-field corrected van der Waals potentials and forces for unpolarized ground-state 
atoms placed within a magnetoelectric medium of arbitrary size and shape. We start from general 
expressions for the van der Waals potentials in terms of the (classical) Green tensor of the electro- 
magnetic field and the atomic polarizability and incorporate the local-field correction by means of 
the real-cavity model. In this context, special emphasis is given to the decomposition of the Green 
tensor into a medium part multiplied by a global local-field correction factor and, in the single-atom 
case, a part that only depends on the cavity characteristics. The result is used to derive general 
formulas for the local- field corrected van der Waals potentials and forces. As an application, we cal- 
culate the van der Waals potential between two ground-state atoms placed within magnetoelectric 
bulk material. 

PACS numbers: 12.20.-m, 42.50.Ct, 34.20.-b, 42.50.Nn, 
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I. INTRODUCTION 

Van der Waals (vdW) forces are non-vanishing forces 
between neutral, unpolarized objects that arise as a con- 
sequence of quantum ground-state fluctuations of the 
electromagnetic field and the interacting objects. The 
first theory of vdW forces within the framework of full 
quantum electrodynamics (QED) goes back to the pio- 
neering work of Casimir and Polder in 1948 Using a 
normal- mode expansion of the quantized electromagnetic 
field, they derived the vdW force on a ground-state atom 
in front of a perfectly conducting plate as well as that 
between two ground-state atoms in free space. 

Whereas the normal-mode technique requires specifica- 
tion of the geometry of the system at a very early point of 
the calculation, alternative approaches based on linear re- 
sponse theory allow for obtaining geometry-independent 
results. Employing the dissipation-fluctuation theorem, 
geometry-independent expressions for the vdW poten- 
tials of one and two ground-state atoms in terms of 
the (classical) Green tensor of the macroscopic Maxwell 
equations and the polarization of the atom can be derived 
[2 i, i, i, i, 0]— results that can be confirmed @, i, El 
within the framework of exact macroscopic QED in lin- 
ear, causal media [Til li~2l. Ii~3j . 

The abovementioned methods apply to atoms in free 
space only, because the local electromagnetic field acting 
on an atom situated in a medium differs from the macro- 
scopic one due to the influence of the local medium en- 
vironment of the atom. Within macroscopic QED, one 
approach to overcome this problem and to account for the 
local-field correction is the real-cavity model [3], where 
the guest atom is assumed to be surrounded by an small, 
empty, spherical cavity. The change of the electric field 
due to the presence of the cavity is then expressed via 
the associated Green tensor that characterizes the local 



environment of the guest atom. 

The real-cavity model has been applied to the spon- 
taneous decay of an excited atom in a bulk medium 
[HI, [HI and at the center of a homogeneous sphere [r| • 
The local-field corrected decay rate for the latter case 
was found to be the uncorrected rate multiplied by a 
global factor and shifted by a constant term. It was 
suggested that the result, reformulated in terms of the 
associated Green tensor, remains valid for arbitrary ge- 
ometries. This was proven correct in Ref. (l71 |. where it it 
is shown that the local-field corrected Green tensor can 
under very general assumptions be decomposed into the 
sum of the uncorrected Green tensor modified by a factor 
and a translationally invariant term. 

Such a decomposition of the Green tensor is not only 
appropriate to obtain the abovementioned results for the 
rate of spontaneous decay. It can also be used to calculate 
local-field corrected vdW potentials. A first hint was 
given in Ref. [Hj], where the local-field corrected vdW 
interaction of two ground-state atoms placed in adjacent 
semi-infinite magnetoelectric media is studied. 

Van der Waals interactions play an important role in 
various fields, where often the interacting objects are em- 
bedded in media so that local-field effects may be ex- 
pected to play a major role. For example, in colloid sci- 
ence the mutual vdW attractions of colloidal particles 
suspended in a liquid influence the stability of such sus- 
pensions [3, El . In particular, clustering of particles 
(flocculation) may occur unless the suspension is suffi- 
ciently balanced by electrostatic repulsive forces 22. 231. 
Examples of vdW interactions of microobjects embedded 
in media can also be found in biology [U [H| , such as 
cell-cell, cell-substratum, cell-virus, and cell-vesicle in- 
teractions [11] , where the investigations are usually based 
on the works of Lifshitz [H, H3] and Hamaker [28| . So, it 
has been studied whether the theory of vdW interactions 
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can explain the phenomenon of biomolecular organiza- 
tion [29[ . In particular, electrostatic and vdW forces 
between a protein crystal and a molecule in solution 
have been (theoretically) compared [30], and the non- 
additivity of vdW interactions between two layers within 
a multilayer assembly such as lipid-water has been in- 
vestigated [3l|. Needless to say that vdW interactions 
represent only one effect among others and that the sys- 
tems are usually very complex, involving a broad range of 
issues such as surface charges, simultaneous existence of 
many different molecules, flexibility of biological macro- 
molecules, fluidity of membranes, and various types of 
non-covalent forces, e.g., vdW, electrostatic, solvation, 
steric, entropic, and structural ones [H, |33|. 

In the present paper, single- and many-atom vdW 
interactions of ground-state atoms embedded in media 
are studied within the framework of macroscopic QED, 
where the local-field correction is included via the real- 
cavity model. The paper is organized as follows. In 
Sec. UU basic formulas for the vdW potentials of ground- 
state atoms are reviewed. Section [TTT] is concerned with 
the real-cavity model, where a general decomposition of 
the local-field corrected Green tensors entering the single- 
and two-atom vdW potentials is presented. The results 
are combined in Sec. IIVI with the general formulas for 
the vdW potential as given in Sec. [ill Closed expressions 
for the local-field corrected single- and many-atom vdW 
potentials are presented and briefly discussed. As an ap- 
plication, the vdW potential is calculated for two atoms 
placed within a magnetoelectric bulk material. Some 
concluding remarks are given in Sec. fVl 



II. BASIC EQUATIONS 

The vdW potential of a (non-magnetic) atom can be 
identified with the position-dependent part of the shift 
of the (unperturbed) eigenenergy, which arises from the 
interaction of the atom with the medium-assisted elec- 
tromagnetic field in the ground state. In leading-order 
perturbation theory, the vdW potential of a ground- 
state atom A at position r A in the presence of a linearly 
responding, dispersing, and absorbing medium, can be 
given in the form 

U(r A ) = ^ |°°duu 2 Tr [ot A (iu)-G^(r A ,T A ,iu) 

(1) 

(for a derivation within the framework of macroscopic 
QED, see, e.g., Ref. Q). Here, 



a A (ui) = lim - 



^ A A a A 



(w A } 2 - <-< 2 
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(2) 



is the atomic ground-state polarizability tensor in low- 
est non-vanishing order of perturbation theory (see, e.g., 
Ref. P), with u k A = {E£ - E$)/h and d| = {i\d A \j), re- 
spectively, being the (unperturbed) atomic transition fre- 
quencies and the atomic electric-dipole transition matrix 



elements. Further, the scattering part of the (classical) 
Green tensor G' 1 ^ (r, r', uj) of the electromagnetic field ac- 
counts for scattering at inhomogeneities of the medium. 
Recall that the Green tensor of the entire system can be 
decomposed according to 



G(r,r',w) 



:(0), 



r,r» + G (1) (r,r» 



(3) 



[G^°\r , r' , oj) , bulk part]. In particular for locally re- 
sponding, isotropic, magnetoelectric matter, the Green 
tensor satisfies the differential equation 



V x 



1 



-V x js{r,uj) 



G(r,r',w) = 5(r-r')l 
(4) 



/z(r,cj) 

(I, unit tensor) together with the boundary condition 

G(r,r',w)— >0 for |r — r'|— >oo. (5) 



The (conservative) force F(v A ) associated with the po- 
tential (fTJ) can be obtained according to 



F(r A ) = -V A U(r A ). 



(6) 



Next, consider two neutral, unpolarized (non-mag- 
netic) ground-state atoms A and B in the presence of 
a linearly responding medium. The force acting on atom 
A due to the presence of atom B can be derived from 
the two-atom potential obtained in fourth-order pertur- 
bation theory as (cf. Ref. [9j) 



U{r A .r B ) 



duu 



2tt jo 

x Tr [a A (iu)-G(r A ,r B ,iu)-at B (iu)-G(r B ,r A ,iu)] . (7) 

Note that in contrast to Eq. (TTJ), the full Green tensor 
G(r A , r B , u) comes into play. The force ¥ A ( B ){r A , r B ) 
on the atom A{B) due to the presence of the atom B(A) 
can be calculated according to 



~F A ( B ){r Al r B ) = -V AiB) U(r Al r B ). 



(8) 



Equation ([7]) can be generalized to the iV-atom potential 
U(r%, . . . ,rAr). In particular, for spherically symmetric 
atoms (ol a = a A [), one finds [35( 



C/(ri, . . . ,rjv) 

(1 + 5 2 n)^ ,/n 



duu ai(iu) • ■ • ajv(iu) 



x STr [G(ri, r 2 , iu) ■ ■ ■ G{r N , n, iu)] , (9) 

where the symbol S introduces symmetrization with re- 
spect to the atomic positions rj., . . . ,rjy, and the corre- 
sponding force on the atom J (J = 1, 2, ... , N) is 



Fj(ri 



,r/,...,rjv) = -VjU(r 1 ,...,r J ,...,r N ). 



(10) 
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III. REAL-CAVITY MODEL 

Let us consider a guest atom in a host medium. In 
this case, the local electromagnetic field at the position 
of the atom, i.e., the field the atom interacts with is dif- 
ferent from the macroscopic field. This difference can 
be taken into account by appropriately correcting the 
macroscopic field to obtain the local field relevant for 
the atom-field interaction. Within macroscopic electro- 
dynamics, a way to introduce local-field corrections is 
offered by the real-cavity model. Under the assump- 
tion that the guest atom (at position r A ) is well sepa- 
rated from the neighboring host atoms, in the real-cavity 
model it is assumed that the guest atom is located at 
the center of a small free-space region, which in the case 
of an isotropic host medium is modeled by a spherical 
cavity of radius R c . Accordingly, the system of the un- 
perturbed host medium of permittivity e(r, ui) and per- 
meability /i(r, lo) is replaced with the system whose per- 
mittivity £ioc(r,^) and permeability /ii oc (r,w) read 

/ n / n J 1 \i \t - r A \< Rc, 

I e(r,L0),fjL(r,L0) if |r - r A \ > Rc- 

(11) 

The cavity radius R c is a model parameter representing 
an average distance from the atom to the nearest neigh- 
boring atoms constituting the host medium; it has to be 
determined from other (preferably microscopic) calcula- 
tions or experiments. Note that on the length scale of 
of the order of magnitude of R c , the unperturbed host 
medium can be regarded as being homogeneous from the 
point of view of macroscopic electrodynamics, i.e., 

e(r,w) = e(r A ,w) \ 
i s i x \ for r-r^t < {l + v )R C , (12) 

with v being a small positive number. 

In the real-cavity model, the local electromagnetic field 
the guest atom interacts with is the field that obeys the 
macroscopic Maxwell equations, with the permittivity 
£i oc (r, oj) and permeability /ii oc (r, to) as given by Eq. (fTTj) . 
Hence, when calculating the vdW potential of the atom 
by means of Eq. |T]), the Green tensor therein is obvi- 
ously the Green tensor of the macroscopic Maxwell equa- 
tions with £i oc (r, ui) and jUi oc (r, cj), respectively, in place 
of e(r, uj) and fi(r, uj), i.e, the Green tensor for the electro- 
magnetic field in the medium disturbed by a cavity-like, 
small free-space inhomogeneity. Clearly, when consider- 
ing the mutual vdW interaction of two or more than two 
guest atoms, each atom has to be thought of as being 
located at the center of a small free-space cavity, and the 
Green tensors in Eqs. © and © are the Green tensors of 
the macroscopic Maxwell equations for the electromag- 
netic field in the medium disturbed by the corresponding 
cavity-like, small free-space inhomogeneities. 

The task now consists in determining such a Green 
tensor. The electromagnetic field inside and outside each 
cavity exactly solves the macroscopic Maxwell equations, 




FIG. 1: Schematic illustration of the interaction of the electro- 
magnetic field with an atom at the center of an empty spher- 
ical cavity (at ta), embedded in a (isotropic) host medium of 
arbitrary size and shape. Simple examples of the three differ- 
ent kinds of processes are sketched. The corresponding Green 
tensors are G^ (ta, va, ui), accounting for all scattering pro- 
cesses within the cavity, G% {r A , ta,w), describing (multiple) 
transmission through the cavity surface and accounting for 
all scattering processes in the host medium excluding reflec- 
tion at the outer boundary of the cavity, and G3 (ta, ta, ui), 
including only the scattering processes involving (multiple) 
reflection at the outer boundary of the cavity. 



together with the standard boundary conditions at the 
surface of the cavity. For an arbitrary current distribu- 
tion j(r, uj), the electric field E(r, t), for instance, can be 
given by means of the Green tensor G(r, r', uj) according 
to the relation 

E(r,i)=j/z J dwwe -1 "' / *dV 'G(r,r' ',w)-j(r' ,w) + c.c. 

(13) 

We begin with the case where a single guest atom is 
embedded in a host medium of finite size and arbitrary 
shape, and, at a first stage, we assume that the host 
medium outside the cavity which contains the guest atom 
is homogeneous — an assumption that will be dropped 
later. According to Eq. ([13]) . the scattering part of the 
Green tensor, G^^r^, r A , uj), describing the electromag- 
netic field reaching the point r A where it has originated 
can be decomposed into three parts, 

G (1) (r A ,i\4,w) = G^\r A ,r A ,uj) + G^ 1 ' (v A , r A , uj) 

+ G£\r A ,r A ,u), (14) 

where Gi(t a ,t a ,u) accounts for (multiple) scattering 

at the inner surface of the cavity, G2 (r A ,r A ,u>) de- 
scribes (multiple) transmission through the surface of 
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the cavity without scattering at the outer surface of the 
cavity, and (r^, r A , oj) is associated with (multiple) 
transmission through the surface of the cavity via scatter- 
ing at the outer surface of the cavity. The decomposition 
is schematically illustrated in Fig. [1] showing a simple 
example each. 

The term (r A , r A , u) in Eq. (fl"4"l) is nothing but 
the scattering part of the Green tensor at the center of 
an empty sphere surrounded by an infinitely extended 
medium. In the case of a magnetoelectric medium of 
permittivity e(u>) and permeability [i(u>), it reads [36j ] 



medium outside the cavity in the case of an infinitely 
extended medium. For a magnetoelectric medium, the 
Green tensor G tr ans(r, w ) reads [36| 

u(u))k 

G t ran S (r,rA, w) = — j D(bj) 



-in 



x [a{q)\ - b(q)vv] e lq (21) 



[k = n{u)cu/c, a(q) = l/q + i/q 2 - l/q 3 , b(q) = l/q + 3i/q 2 
- 3/g 3 , q = [n(io)u>/c]\r - r A \, v = (r - r A )/(\r - r A |)], 
where 



G[ 1) (r A ,r A ,u J ) = £-C(u;)\, 



(15) 





zh^\z) 


-e{u)h^(z) 


z h^\z ) 


e(uj)h[ 1 \z) [z ji(zo)]' - ji(z ) 


zh { ?\z) 


i 



where 



C(u) = 



(16) 

[zq = u)R c /c, z — n(uj)zQ, n(ui) — ^Je(uj)h(uj) , the primes 

indicate derivatives], with j\(x) and hi (x), respectively, 
being the first spherical Bessel function and the first 
spherical Hankcl function of the first kind, 



. . sin(x) cos(x) 
3i{x) = — j 



(17) 
(18) 



At this point it should already be mentioned that since 
the vdW potential depends on the Green tensor at all 
frequencies, an expansion for small \lu\R c /c of the form 



CM 



e{ui) - 1 c 3 



2e(uj) + 1 iu z Rl 
9 e 2 (u))[5fj,(ui) - 1] -3e(w) - 1 
5 
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[2e(w) + l] 2 
e(a;)n 3 (a;) 

[27RTTp 



-l + O 



iu}R c 

UjRr 



(19) 



and neglect of the term 0(ujR c /c), as it is done in 
Ref. [17| for the calculation of the decay rate of an excited 
atom, has to be applied with care (see Sec. IIV A2p . 
Following the arguments given in Ref. [17| . the term 

G2 \r A ,r A ,u) in Eq. (fT4|) can always be written in the 
form 



G$\t a ,t a , U ) = / 2 HG^(r A ,r,, W ), (20) 

where G[^ cd (r A ,r A , lu) is the scattering part of the Green 
tensor of the undisturbed host medium (i.e., the medium 
without the cavity), and the global factor f(u>) can sim- 
ply be determined from examining the Green tensor 
Gtrans( r ; r A, which links a source at the center of the 
cavity with the transmitted field at a point r in the 



(i) 



D(w) = 



JiOo) 


z hY >(zo) 


1 (1) 

- [zoji(za)}' h\ >(z ) 


mH{j'i(2o) 


zh [ i\z) 


- [z ji{z )]' \z)\ 



(22) 

Comparing Eq. (j2"Tj) with the Green tensor for an undis- 
turbed, infinitely extended host medium, 

G^(r,r x , W ) = ^ [a(q)\ - %)vv] a* (23) 

we see that the factor f(ui) is equal to D(ui), 

f(u) = D(u). (24) 
Inspection of Eq. (|22|) for small \ui\R c /c shows that 



D(u) 



3e{u) 
2e(w) + 1 



O 



UjRr 



(25) 



Note that in leading order, the factor D(cu) is the same 
as that found for a purely electric system (cf. Ref. [ijj )• 
As shown in Ref. [l7|, the term G^ (r A ,r A ,u) in 
Eq. (jT4j) behaves like (luR c /c) 3 for a sufficiently small 
cavity in a dielectric medium. The same asymptotic be- 
havior is observed for magnetoelectric media, since, for 
sufficiently small R c , the leading-order contribution to 
G3 (t a ,t a ,u>) is determined by the electric properties. 
Recalling that in the real-cavity model the cavity radius 
is assumed to be small compared with the relevant atomic 
and medium wavelengths, G^\r A ,r A ,uj) is in any case 

small compared with G± (r A ,r A ,uj) and G^\r A ,r A ,uj) 
and can thus be neglected, so that Eq. (fl"4|) [together with 
Eqs. ([IS}. (f2T)]) . and (|2"I))] takes the form 



G W (r A ,r A ,Lu) 



tul 

Gtt 



C( W )l+^HG« d (r A ,r A , W ) 



(26) 

So far we have concentrated on homogeneous host me- 
dia of arbitrary shapes. An inhomogeneous medium can 
always be divided into a (small) homogeneous part that 
contains the cavity, such that the condition (fl"2)) is satis- 
fied, and an inhomogeneous rest. Since all the scattering 
processes resulting from the inhomogeneous rest can be 
thought of as being taken into account by the scattering 



l r ) 




FIG. 2: Schematic illustration of the interaction of the electro- 
magnetic field with two atoms, each placed in a spherical cav- 
ity, surrounded by a host medium of arbitrary size and shape. 
The relevant Green tensors G4(rA,rA,^) and Gs(rA,rA,^) 
account for (multiple) transmission through the cavity sur- 
faces and (multiple) scattering at the inhomogeneities of the 
host medium including (multiple) scattering at the outer 
boundaries of the cavities, respectively. 



sion through the surfaces of the cavities such that the 
field at r outside the cavities which results from yb, 
may be calculated from Ds(w)G me d(r, ?b, to). Due to 
the property G(r, r',w) = G (r',r, a;) of the Green ten- 
sor, we also find that Z?A(^)G mc d(r J 4, r, oj) determines 
the electromagnetic field at point ya originating at point 
r outside the cavities. We may thus write G(r a, yb,w) 
in the form 

G(r A ,r B ,u) = D A {uj)G mcd (YA,Y B ,uj)DB{uj). (28) 



IV. LOCAL-FIELD CORRECTED VAN DER 
WAALS INTERACTION 

A. Single atom within an arbitrary 
magnetoelectric medium 

Within the framework of the real-cavity model, local- 
field corrected vdW interactions can now be calculated by 
combining the formulas given in Sec.HTland lllli where we 
again begin with the case of a single guest atom. We sub- 
stitute G^(ta,ya,lj) as given by Eq. (|26j) into Eq. ([T|) 
to obtain 

U(t a ) = U 1 (t a ) + U 2 (t a ), (29) 

where 



part of the Green tensor of the host medium, Eq. (|26|) 
still holds in this more general case. Note that a homo- 
geneous medium of finite size can already be regarded as 
a special case of an inhomogeneous system. Hence, e(ui) 
and fJ>(u>), respectively, can be safely replaced by e a (lj) 
= e(r A ,u)) and /j,a{w) = /u(r^,a;) in the equations given 
above [C(u>) i— > Ca{u), D(ui) h- * D a (cj)]. 

Next we consider two atoms within a host medium of 
arbitrary size and shape, both located at the center of 
a small cavity, such that Eqs. (jTTJ) and (fT2|) hold. In 
this case, there are two contributions to the Green tensor 
G(ya,yb, oj) that determines the electromagnetic field at 
the position r A of atom A, originating at where is 
atom B located, 

G(i\4,r B ,w) = G±{va,t b ,oS) + G 5 {y a ,y b ,oj). (27) 

Here, G^yatYbjOj) accounts for (multiple) transmission 
through the cavity surfaces and (multiple) scattering at 
the inhomogeneities of the medium, excluding scatter- 
ing at the outer boundaries of the cavities. As in the 
single-atom case, processes that involve (multiple) scat- 
tering at the outer surfaces of the cavities, described by 
G$(ya, yb, oj), are neglected. The situation is sketched in 
Fig.H 

By means of the same arguments as in the single-atom 
case, we find a global factor that accounts for transmis- 



Ux{y a ) 
U 2 {v a ) 







12tt 2 c 



duu 3 CA{iu)Tr a A (iu), (30) 



2tt Jo 



duu D A (iu) 



x Tr 



a A {iu)-G { l d (rA,r A ,iu) 



(31) 



[C A (iu) and Da(iu) according to Eqs. (fT?J)) and (|2"2"j) . re- 
spectively] . 



1. Linear approximation 

In order to make contact with a microscopic descrip- 
tion and provide a deeper insight into the validity of the 
real-cavity model, let us first restrict our attention to 
to weakly polarizable and/or magnetizable host material 
so that the linear approximation applies. In this case, 
we may regard e(r,a>) and fj,(r,u>) as being only slightly 
perturbed from their free-space values, 



e(r, u) 



•X(r,w), 
■C(r,w) 



(32) 
(33) 
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[|x(r,w)| < 1, \((r,oj)\ < 1], and approximate Ui(r A ), 
Eq. ([!)]), by 



with 



Ui(r A ) = --^r c J duu 3 a A (i 



-2uR c /c 



2c 2 



(34) 



where we have assumed a spherically symmetric atom. 
By straightforward calculation, one can then show that 

U 1 {v A )=p{r A ) f d 3 sU{v A ,s), (35) 

J\ S -r A \>R c 

where p(r) is the number density of the medium atoms, 

U(r A , s) = U c \r A} s) + f/ mag (r A , s) (36) 

(r^^s) is the vdW potential of two unpolarized, neutral 
(ground state) atoms in free space, with the electric and 
magnetic parts being given by [T(| 

ha 2 f°° 

U (r A ,r) = - / duu A a A {iu)a(iu) 

2vr J 

iu 

> 

dua A {iu)a{iu)g(ul/c) (37) 



x Tr[Gf ree (r,4, r, iu)- G licc {r, r A , iu)} 
h 



327T 3 £2/6 Ju 



and 



L7™*(r A ,r)=-M HAH 
27T Jo 



u 2 a A (iu)P(iu) 



Tr| G free (rA,s,iii)xV : 



•V s xG 



free 



(s,r^,i«) \ 



32tt 3 Z 4 



duu a A (iu)(3(iu)h(ul/c) (38) 



\g(x) = 2e- 2x {3 + 6x + 5x 2 + 2x 3 +x 4 ), h{x) = 2 e - 2:E (l + 
2a; + a; 2 ), ^ = |r,4 — r|; Gf ree (r, s, iu), free-space Green 
tensor; a(iu), polarizability; f3(iu), magnetizability of the 
medium atoms] , where we have made the identifications 



X(r,w) = e 1 p(r)a(u), 
C(r,w) = n Q p(r)/3(iu). 



(39) 
(40) 



To find the linear approximation to U2(r A ), Eq. (f3"Tj) . 
let us first consider 

G mtd( r A, T A,u) = G med (r A ,r A ,u) - G { ^ ed (r A ,r A ,w). 

(41) 

Applying the linear Born expansion (see, e.g., Ref. [l7| ) 

G(r, r', w) = G froc (r, r', w) + AG(r, r', u) (42) 
AG(r, r', u) = AG cl (r, r', to) + AG mag (r, r', w), (43) 



AG ol (r,r» = ^ / d 3 sx(s,^) 

x G froc (r,s, w)-Gf ree (s,r',w), (44) 

AG mag (r,r» = - / d 3 sC(s,w) 

x G froc (r,s,w)x V s •V s xGf re e(s,r',w) (45) 



to the calculation of the Green functions G mc d(r.4, Fa, w) 
;(0) / 



and G^ d (r j4 ,r j4 , w) in Eq. (gT|), we derive 



G mod( r ^> r ^^) 



c ./V 



/ d 3 s x(s, w)G bee (r A , s, w)-G free (s, r^, w) 



d 3 sC(s,w) G frcc (r,s,w)xV s ■V s xG fa (s,r',w) 

w 2 x(r j4 ,w) f ,3 _ , v _ , , 
-2 / d J sG fr ee(rA,s,w)-Gf rco (s,r j4 ,c i ;) 



+ C(rA,^) / d 3 s G frcc (r,s,w)xV s V s xGf Iee (s, r', 10 



(46) 

(V, volume of the unperturbed host medium). Equation 
(|46|) reveals that, in contrast to D a (oj), G^l d (r A ,r A ,Lu) 
has no zero-order contribution with respect to x( r > w ) 
and C(r, w). Thus, the linear approximation to U2(j A ) 
can be obtained by substituting Eq. ([4*6")) into Eq. (pTTj) . 
letting I?a(w) = 1, and applying Eqs. (j39j) and (|40|) . which 
results in 

tMrA) = / d 3 sp(s)[/(r j4 ,s) 
iv 

-p(r A ) Jd 3 sU(r A ,s). (47) 

As expected, the overall potential in linear approximation 
is simply 



U(v A ) = U 1 (r A ) + U 2 (r A ) 

= [ d 3 S p(s)(7(r A ,s) 
Jv 



(48) 



(V , volume of the host medium excluding the volume of 
the cavity), where we have recalled Eq. (fT2")l . 

The results found in linear approximation, particularly 
Eq. (|48p. indicate that R c can indeed be identified with 
an average inter-atomic distance. Hence, the real-cavity 
model can be regarded as being valid down to a micro- 
scopic scale. Since ^(r^) can be described as the sum of 
the vdW interactions between the guest atom and all the 
atoms in an infinitely extended homogeneous medium, it 
does not lead to a force but to a more or less significant 
shift of the overall potential, which is dominated by the 
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interaction of the guest atom with the neighboring host 
atoms. Note that we have demonstrated these results 
up to linear order in x(r,uj) and £(r, ui). In principle, a 
generalization beyond the linear Born expansion should 
be possible, although associated with extensive calcula- 
tions. In this case many-atom vdW interactions have to 
be taken into account. 



2. 



ond the linear approximation 



Let us return to the integrals in Eqs. (|30|) and (|3Tj) . To 
find the local-field correction, one has in fact to consider 
the leading-order terms in R c after the frequency inte- 
grals have been performed. However, since e A (u>) — 1, 
fj, A (u>) — 1, and a A (uj) provide cut-off functions for the it- 
integrals [giving their main contribution for u < w max , 
with w max being the maximum of all relevant atomic 
and medium resonance frequencies] and assuming that 
Wmax-Rc/c < 1, we may expand C A (iu) and D A (iu) in 
the integrands, keep the leading-order terms in uR c /c, 
and integrate over u afterwards. This procedure can be 
justified by the good convergence behavior in the asymp- 
totic limit of small R c , which we have confirmed numer- 
ically. 

Thus, using Eq. (19]), we can write Ui(r A ), Eq. (|30)l . 
in the form 



U^ta) 



du 



Sa(iu) — 1 c 2 
2e A (iu) + 1 ~Rj 



s\(iu) [1 — 5fi A (iu)] + 3eA(iu) + 1 u 2 



5[2e A (iu) + lY 



xTr a A (iu). (49) 



As already mentioned, Ui(r A ) can be ignored if one is 
only interested in the force on the atom (see also App.JA}. 
Similarly, using Eq. (|25|) . we can write U^ta), Eq. (f3Tj) . 
as 



U 2 (r A ) = 



2ir 



duu 



3eA(iu) 



n 2 



2e A (iu) + 1 



x Tr a A («it)-G^L(r A ,r A ,w) . (50) 



With Eq. © and the condition (fT2|) . one obtains 

-i 2 



F(r^) 



2tt 



dtttt 



3£a(»m) 



2£a(« 



1 



x V A Tr q a (jm) • G { ^l d (r A ,r A ,iu) (51) 



for the local-field corrected vdW force on a ground-state 
atom within an arbitrary magnetoelectric system. 

In the limit of the atom being embedded in a very di- 
lute medium, e A (iu) -1<1 and u A (iu) -1<1, the 
local-field correction becomes negligible and Eqs. (|49p- 
(ISTj) approach the uncorrected result [Eqs. ((T|) and ©], 



as expected. As a note of caution, we remark that, quite 
contrary, in the limit R c — -> 0, the local-field corrected re- 
sult does not approach the uncorrected one, which might 
be obtained when the atom is thought of as not being 
placed in an empty free-space region, but being placed in 
a region which is continuously filled with medium. This 
is due to the fact that there is no continuous R c — > tran- 
sition from the corrected setup (an atom inside an empty 
free-space cavity surrounded by a medium) and the un- 
corrected one (an atom inside the medium): Even as the 
cavity radius is made smaller and smaller, the atom will 
always remain inside a free-space region. 

Equation (|5ip shows that, in comparison to the uncor- 
rected vdW force ([6]) , the contributions to the it-integral 
are enhanced by a factor between 1 and 9/4 owing to 
the local-field correction. Note that the enhancement 
factor varies with u and the uncorrected contributions 
can in general (in particular, for magnetoelectric media) 
have different signs for different values of u. Hence, the 
local-field corrected vdW force can be enhanced or re- 
duced. For example, for (purely electric) media impor- 
tant in biological applications, local-field corrections can 
lead to enhancements by factors of 2.22 [for water with 
e(Q) ss 80] or 1.86 . . . 2.19 (for typical plasma membrane 
compositions with e(0) w 5 ... 40 |24|). 



B. Two or more atoms within an arbitrary 
magnetoelectric medium 

Using Eqs. ^ and (|28p. the local-field corrected two- 
atom vdW potential reads 



U(r A ,r B ) 



2tt 



duu Da (iu)Dg (iu) 



Ty[ola (iu) ■ Gmcd (r a ,r B ,iu)-a B (iu) ■ G mo d (re , r A , iu)]. 

(52) 



In analogy to Sec. lIV A~2l we may expand the coefficients 
D A (iu) and Db(«m) to leading order in uR c /c, as given 
by Eq. (f2"5| . and integrate afterwards. Hence, Eq. (|52"|) 
leads to 



U(r A ,r B ) = -— 



/>oo 

x / duu 4 
Jo 



3e A (iu) 
2e A (i 



3eb(iu) 



-1 2 



2e B (i 



xTr[a A (iu)-G m ^{T A ,TB,iu)-a B {iu)-G mB i(TB,r A ,iu)]. 

(53) 

In comparison to the uncorrected potential (|7|) , the con- 
tributions to the it-integral are enhanced by a factor 
which varies between 1 and 81/16, so that altogether an 
enhancement or an reduction of the local-field corrected 
two-atom force may be observed (recall the remark at the 
end of Sec. lIV A 211 . 
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Equation (|53[) can easily be extended to the mutual 
interaction of N atoms. Substituting Eqs. (j2"5)) and l|28p 
into Eq. J9}, one finds that, for isotropic atoms, the local- 
field corrected iV-atom potential is given by 



(1 + S 2 n)tt 



duu 2N 



3ei(m) 


2 


3ejv(iw,) 


_2e 1 {iu) + 1 




2Sfq(iu) + 1 



c*i(iu) ■ ■ ■ ajv(iw) 
x Tr[G m ed(ri,r 2 ,m)---G ro ed(rAr,ri,m)]. (54) 



C. Example: Two atoms within a bulk 
magnetoelectric medium 

As a simple example, let us consider two isotropic 
ground-state atoms embedded in magnetoelectric bulk 
material characterized by e(to) and n(u>). Substituting 
Eq. (J53J into Eq. (|53"]l. we obtain 



U(r A ,r B ) = - 
3e(«u) 



ft 



aA(iu)a.B{iu) 
e 2 (iu) 



{3 + 6n(iu)ul/c + 5{n(iu)ul/c 



2e(iu) + 1 

f 2[n{iu)ul/cf + [n(iu)ul/c] 4 } e ~^(™)ui/c^ (55) 



where I — \r A — rg|. In the retarded limit, I 3> c/w m i n 
(with w m in denoting the minimum of all the relevant 
atomic and medium resonance frequencies), the main 
contributions to the u-integral come from the range 
where 



aA(B){iu) ~ au(s)(0), e(iu) ~ e(0), n(iu) ~ //(O) 
With this inserted, Eq. (|55|) simplifies to 

U{r A ,r B ) = — ^, 



where 



23ftc QA(0)aB(0) 
647r 3 £ 2 n(0)e 2 (0) 



3e(0) 



1 4 



2e(0) + 1 



(56) 
(57) 

(58) 



In the non-retarded limit, !<c/n(0)w ma st) we may set 
g -2n(iu)ul/c ^ i anc j approximate the term in the curly 
brackets by 3. In this way, Eq. (|55|) reduces to 



U(ta,t b ) = - 



Cnr 
16" 



(59) 



where 

Cnr 



3ft 



-o Jo 



00 aA(iu)aB(iu) 
au 



s 2 (iu) 



3e(iu) 



2e(iu) + 1 



(60) 

Equations (|F7|) and ([55^) generalize the well-known re- 
sults for the uncorrected mutual vdW interaction of two 



atoms in bulk material [2jJ. We see that the local-field 
corrected potential follows the usual l~ 7 and l~ 6 depen- 
dence in the retarded and non-retarded limits, respec- 
tively. However, due to the influence of the local-field cor- 
rection, the potential/force can be noticeably enhanced 
in dense media. 



V. SUMMARY 

In this paper, we have studied the local-field corrected 
vdW interaction of ground-state atoms embedded in a fi- 
nite magnetoelectric medium of arbitrary size and shape. 
We have started from the familiar expressions for the one- 
and many-atom vdW potentials in leading-order pertur- 
bation theory, given in terms of the associated (classical) 
Green tensors and the atomic polarizabilities. Employ- 
ing the real-cavity model for the local-field correction and 
using the methods developed in Ref. [I3|, we have de- 
rived general expressions for the required local-field cor- 
rected Green tensors. Using these results, we have ob- 
tained geometry-independent, general formulas for the 
local-field corrected one- and many-atom (ground-state) 
vdW potentials and the corresponding forces. 

The derivation is based on the assumption that the 
radius of the free-space cavity, in which a guest atom is 
thought of as being located and which is a measure of the 
average distance of the guest atom from the surrounding 
host atoms, is small compared to the other characteristic 
lengths of the problem. That is to say, it should be small 
compared to the distance of a guest atom from surfaces of 
discontinuity of the host medium, the distances between 
the guest atoms in the case of many- atom interactions, 
and the wavelengths associated with the electromagnetic 
response of the guest atoms and the host medium. It 
should be pointed out that under these conditions the 
forces do not explicitly depend on the cavity radius. 

We have found that the local-field correction leads to 
a position-independent shift of the single-atom potential, 
while at the same time enhancing the contributions to 
the position-dependent part by a frequency-dependent 
factor that may vary between 1 and 9/4. The frequency- 
dependent contributions to the two-atom potential, on 
the other hand, are enhanced by a factor between 1 and 
81/16 owing to the local-field correction. The contribu- 
tions to the general iV-atom potential are enhanced ac- 
cordingly. In all cases, the respective vdW force can be 
enhanced or reduced. 

As a simple example, we have calculated the mutual 
vdW potential of two ground-state atoms placed within 
a magnetoelectric bulk medium. The result shows that 
the potential follows the well-known asymptotic l~ 7 and 
l~ 6 power laws in the retarded and non-retarded regimes, 
respectively, where the local-field correction leads to an 
enhancement of the associated proportionality constants. 

We conclude by noting that, besides vdW interaction 
and spontaneous decay, the local-field corrected Green 
tensors obtained can also be used to account for the ef- 
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fects of the local field in a number of important processes 
(e.g., fluorescence, light scattering, energy transfer, etc.). 
Evidently, one can expect that these processes are af- 
fected by local-field corrections in a similar way as found 
here. 
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APPENDIX A: ATOM INSIDE A FREE-SPACE 
CAVITY 



we have 



G ( 1 )(r A ,r^ W ) = £^ + l)Cf( 1 



1=1 



A-l 
$A 



(21 + 1)1! 



x [l(l + l)e r e r + \(l + lf(e e eg + e e )] (A3) 

for the small £4 expansion of the Green tensor. As seen, 
in the limit £4 — > 0, only the I = 1 term survives, leading 
to the potential (|30|) . The force on an atom close to 
the cavity center is determined by the derivative of the 
Green function with respect to ta- Therefore, for small 
displacements of the atom from the cavity center, it is 
given by the I = 2 term in the above expansion, and we 
find the linear force 



In Sec. IIV1 we have encountered the potential term 
Ui(r A ), Eq. (|30[) . which is equal to the potential of an 
atom at the center of a free-space cavity embedded in 
a bulk medium of permittivity ea(>jj) and permeability 
/j,a(w)- More generally, the potential of an atom at a 
small displacement from the center of the cavity is 
given by Eq. (fTJ) with the Green tensor 36] 



G^(r A ,T A)U ) = ^£(2Z + l){cf HZ(Z + 1) 



Dm 



1=1 



+ <TMif(60 



x (e g e g + e^e^Y (Al) 



(e r , eg, e^, spherical unit vectors associated with r^; 
£a = uJTa/c). Here Cf is given by Eq. ([To]) , with the 
subscript 1 — ► /, and Cf 1 = CF[e <-> fi\. From 



(21 + 1)1) 



for £4 < 1, 



(A2) 



with 



F(r A ) = -e r U(r A ) = Kv A (A4) 

OTA 



K = — 



12^ 2 e c 5 7o 



duu 5 a A (iu)C2(iu), (A5) 



[recall Eq. (fl6|) ]. which vanishes at the cavity center. For 
small cavity radii Eq. (|A5j) reduces to 



K = 



127T 2 £ c 5 



duu 5 otA(iu) < 90 



e(iu) — 1 c 5 



75 e 2 (iu)[7n(iu) + 3] - 6e(iu) 
Y 



[3e(iu) + 2] 2 



3s(iu) + 2 u 5 R5 



(A6) 



u 3 m 



leading to a force pointing away from the center of the 
cavity (corresponding to an unstable equilibrium) for ma- 
terials with purely dielectric properties and to a harmonic 
force for purely magnetic materials. 
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